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Abstract 

The two-dimensional Navier-Stokes equations are rewritten as a system of 
coupled nonlinear ordinary differential equations. These equations describe 
the evolution of the moments of an expansion of the vorticity with respect to 
Hermite functions and of the centers of vorticity concentrations. We prove the 
convergence of this expansion and show that in the zero viscosity and zero core 
size limit we formally recover the Helmholtz-Kirchhoff model for the evolution 
of point- vortices. The present expansion systematically incorporates the effects 
of both viscosity and finite vortex core size. We also show that a low-order 
truncation of our expansion leads to the representation of the flow as a system 
of interacting Gaussian (i.e. Oseen) vortices which previous experimental work 
has shown to be an accurate approximation to many important physical flows 

®- 

1 Introduction 

In this paper we represent solutions of the two-dimensional Navier-Stokes equations 
as a system of interacting vortices. This expansion, which generalizes the Helmholtz- 
Kirchhoff model of interacting point vortices in an inviscid fluid, systematically in- 
corporates the effects of both vorticity and finite vortex core size. Furthermore, we 
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give conditions which guarantee the convergence of our expansion. Incompressible 
viscous flow has two standard analytic representations: a formulation in terms of the 
primitive velocity and pressure variables, and a formulation in terms of the veloc- 
ity and vorticity variables |7]. The velocity- vorticity representation has particular 
advantages when boundaries are unimportant, since vorticity cannot be created or 
destroyed in the interior of a fluid. The vorticity field can also be directly related to 
physically observed flow structures such as line and ring vortices. 

In two space dimensions, the vorticity field has the additional advantage of re- 
ducing to a scalar. An early representation of two-dimensional flow in terms of 
moving point vortices was developed by Helmholtz-Kirchhoff [I] and by Helmholtz 
[15] . The point vortex model has been studied extensively - a thorough review of 
the model and recent developments are described in the monograph by Newton [TT] . 
While the Helmholtz-Kirchhoff point vortex model captures many of the basic phys- 
ical phenomena observed in two-dimensional rotational flows, experiments with even 
relatively simple vortex configurations exhibit complications far beyond the point 
vortex predictions [9]. Additionally, the classical point vortex model neglects the 
effects of viscosity. However, these experiments also reinforce the idea that in many 
circumstances the fluid flow may be well approximated by a collection of interacting 
vortices - albeit vortices with finite core size, subject to the effects of viscosity. A few 
recent studies of the interaction of viscous vortices can be found in [lj [5j El [HI [91 [H] . 
The main focus of these papers is the merger of two like signed vortices. In [8] the 
authors use a spatial moment model for 2-D Euler equations which later incorpo- 
rates weak Newtonian viscosity to derive equations of motions for two like signed 
vortices. A metastable state is found before merger which consists of two rotating, 
near-circular, vortices. More recently, it has been conjectured that in a two-vortex 
system the profiles relax to a pair of gaussian vortices before merging [6] . Thus, it is 
of interest to extend and generalize the Helmholtz-Kirchhoff point vortex model to a 
model that incorporates non-zero vortex core size and viscous effects while retaining 
its basic form. Such an extension is the goal of this paper. 

The governing equations for the velocity (u) and pressure (p) variables are 



where p is the fluid density and v is the kinematic viscosity. Taking the curl of ([T}) 
and (T5]) gives 
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w = V x u, V • u} — 0, 



(4) 



which are the governing equations for the velocity-vorticity variables. For two- 
dimensional flows, the vorticity vector uj is perpendicular to the plane of the flow, 
and the third term on the left-hand side of ([3]) vanishes. The condition V • uj = is 
identically satisfied, and equation ([3]) then reduces to the single scalar equation: 



where lo is the single, non-zero component of the vorticity. A drawback of the 
formulation in §5§ is that the velocity of the fluid is still present in the equation. 
However, assuming that the vorticity field is sufficiently localized, the velocity vector 
can be computed in terms of the vorticity uo by the Biot-Savart law 



where for a two-vector z — (zi, z 2 ), z 1 - = (—z 2 , Zi). 

In this paper we use equations fl^D and (jSJ) to develop a vorticity representation 
of two-dimensional viscous flow. Our representation is based on a decomposition of 
the vorticity field into a set of moving distributed vortices. Differential equations 
are derived for the motion of the vortex centers and for the time evolution of the 
vortex distributions. The evolution of each individual vortex is represented as an 
expansion with respect to a sequence of Hermite functions. Such expansions have 
proven useful in theoretical studies of two-dimensional fluid flows ([2], [3]) and the 
leading order term in this expansion is precisely the Gaussian vortex (i.e. Oseen 
vortex [13], [12]) whose utility as an approximation for vortex interaction was shown 
in [9]. We show that the coefficients in this expansion satisfy a system of ordinary 
differential equations whose coefficients can be explicitly represented in terms of a 
fixed, computable kernel function. We also prove the convergence of this expansion. 
It is shown that our representation reduces in the appropriate limit to the Helmholtz- 
Kirchhoff model, and allows at the same time arbitrarily complex evolution and 
interaction of the moving vortices. 

In the present paper we concentrate on the mathematical formulation of the gen- 
eralized Helmholtz-Kirchhoff model. In future work we will explore the predictions 
of this model both numerically and analytically in a number of different physical 
settings. 




(5) 
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2 The "multi-vortex" expansion 

In this section we separate the solution of the vorticity equation into N components 
and derive separate evolution equations for each component. From a physical point of 
view this decomposition will be most useful when each of the components corresponds 
to a localized region of vorticity (e.g. a vortex) well separated from the other lumps 
but the mathematical development described below is well defined without regard to 
these physical considerations. However, with this application in mind, we will often 
refer to each of the components "vortex" . 

Consider the initial value problem for the two-dimensional vorticity equation 



u(x, 0) = ljo(x) 

where u is the velocity field associated to the vorticity field u. We begin by decom- 
posing the initial vorticity distribution by writing 



Of course this decomposition is not unique - even the number of pieces, N, into which 
we decompose the vorticity is up to us to choose. In general the choice we make will 
be motivated by physical considerations, however, for the development below, all we 
require of the decomposition is that the total vorticity of each vortex is non-zero, i.e. 



— — = uAuj — u • Vw, 
dt 

u = uj(x,t), x e M 2 , t > 



(7) 







If (J9]) is satisfied we define Xq by 




(10) 



or equivalent ly 




(11) 



We now write the vorticity for t > as 



N 




(12) 
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and the velocity field as 



N 

u(x,t) = ^2u j (x-x j (t);t) (13) 

3=1 

where u^(y,t) is the velocity field associated to u^(y,t) by the Biot-Savart Law. Of 
course we still have to define the equations of motion for uj^{y,t) and x^(t). 

The centers of the vorticity regions, x^(t), and the vorticity regions themselves 
evolve via a coupled system of ordinary-partial differential equations constructed so 
that in the limit of zero viscosity and when the different components of the vorticity 
happen to be point vortices (i.e. Dirac-delta functions) we recover the Helmholtz- 
Kirchhoff point vortex equations. If we take the partial derivative of (fl2l) and use 
the equation satisfied by the vorticity, we find: 

N N 

d t u(x,t) = ^<W(s - s J '(t),^) • Vu j (x -x j {t),t) (14) 

3=1 3=1 
N N / N \ 

= ^z/Ai^x-z 3 ' (*),*) - ( 5^ — ar € (t), *) J • Wu j {x-x j (t),t) . 

3=1 3=1 \t=l J 

Given this equation it is natural to define as the solution of the equation: 

o j / N \ 

'x-x j (t),t) = vAcu j (x-x j (t),t)- J]u £ (x-x £ (t),t) • W(x-x j (t),t) 



dt 



=1 

+x j (t) ■ V(J{x - x j (t),t) , j = 1, . . . , N. (15) 



To close this system of equations we must specify how the centers of vorticity 
xi(t) evolve. We impose the condition that the first moment of each vorticity region 
must vanish at every time t > 0, i.e. we require that 

(x - x j (t))cu j (x - x j (t),t)dx = for all t > 0, j — 1, . . . N . (16) 

(Note that this equation really contains two conditions - one for each component of 
(x— x e (t)).) We impose this condition to fix the evolution of x j (t) because Gallay and 
Wayne have recently shown [2] that if one considers the evolution of general solutions 
of (J7J) the solution will approach an Oseen vortex, and the rate of the approach will 
be faster if the vorticity distribution has first moment equal to zero. Solutions of 
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(E3) preserve the first moment, and hence if the initial conditions have first moment 
equal to zero the solution will have first moment zero for all time. The equations 
(11511 no longer preserve the first moment and thus we impose this condition for all 
time, which then defines the motion of the center of vorticity. 

Note that if we change variables in f|T6|) to z = x — x-*(t), we find 

zu j {z,t)dz = . (17) 

Since this equation holds for alH > we can differentiate both sides with respect to 
t to obtain 

/ zd t u j {z,t)dz = . (18) 

JR 2 

Using ([TBI) we can insert the formula for d t u>^ into this integral and we obtain: 

uf zAu j (z,t)dz- f z ( YV(z + -x l (t),t) ) ■ Vu 3 {z,t)dz 
Jr 2 Jr 2 J 

+ [ z (x j (t) ■ Vu j {z,t)) dz = . (19) 
Jr 2 

We first note that if we integrate twice by parts, we have 

/ zAu j {z,t)dz = . (20) 
Jr 2 

Next if we take the n th component and integrate by parts we find 

z (x{(t) ■ Vuj j (z, t)) dz = xl(t) [ zd Zn co j {z, t)dz = -rrijiUt) , (21) 



where rrij = f^u^z, t)dz and n — 1, 2. 

Remark 2.1. Note that the equations f|T5|) do preserve the total integral ("mass") 
of the solution so this definition of nij is consistent with (Q . 

Finally, recalling that the velocity field is incompressible we can rewrite 

N \ / N 

J2 ui ( z + ^'W - At),t) ■ vt^'Oz, t) = v ■ lj2 u ^ z + ^'W - x e {t),t)u j {z, t) 



6 



so, again considering the n th component and integrating by parts we have 
I z n [ ^u\z + x j {t) -x e (t),t) ) -Vu j {z,t)dz 

J** \ t=1 J 

= -[ (j2<^ +xj ^- 




Thus, if we combine ( 1201) (l21j) . and (1221 we see that (1X01) reduces to the system 
of ordinary differential equations for the centers of the vorticity distributions: 



dt 



(i) = -LW (j£( z + J(t)-a*(t),t)ui(z,t))dz, (23) 

m j jr[ Jr 2 



supplemented by the initial conditions (jTUl) . while the N components of the vorticity 
evolve according to the partial differential equations (1151) with initial conditions 

u j (z,0) =ui(z + x 3 ) (24) 

obtained by combining (jSJ) and 012|) . 

Remark 2.2. Consider (123]) in the limit in which the components are all point 
vortices, i.e. u>^(z,t) = rrijS(z), with S(z) the Dirac- delta function. Recall that the 
velocity field associated with such a point vortex is 



2 1 

U n (Zi, Z2,t) = — y € n jZj , 2 2 ~r 



where e mj - is the antisymmetric tensor with two indices. Then if we ignore the 
(singular) term with i = j in the sum on the right hand side of (|23p we find 

"*~ = JE, mt — m^w — (25) 

These of course are just Helmholtz-Kirchhoff equations for the inviscid motion of a 
system of point vortices. Thus, our expansion can be regarded as a generalization 
of the this approximation which allows for both nonzero viscosity and vortices of fi- 
nite size. To justify omitting the term with i = j on the right hand side of (j25|) 
we note that if we approximate the delta function with a narrow, Gaussian vortic- 
ity distribution, and u- 7 by the corresponding velocity field, this term will vanish by 
symmetry. 
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3 The moment expansion; case of a single center 

In this section we introduce another idea - an expansion of the vorticity in terms of 
Hermite functions. Then, in the next section we will combine the Hermite expansion 
with the multi-vortex expansion of the previous section. 

The moment expansion is an expansion of the solution of the vorticity equation 
in terms of Hermite functions. Define 

<Mx,t;A) = 4^ |2/A2 ( 26 ) 
where A 2 = A 2 , + 4z4. Three simple facts that we will use repeatedly are 

(i) d t (poo = ^A(/>oo 

(ii) f R2 (j)oo(x,t] X)dx = 1 for all t > . 

(iii) Finally, and crucially for what follows, the vorticity function u(x, t) = a<froo{x, t) 
is an exact solution (called the Oseen, or Lamb, vortex) of the two dimensional 
vorticity equation for all values of a. 

Note that we will often supress the dependence of O o on A when there is no fear of 
confusion. 

We now define the Hermite functions of order (ki, k^) by 

fcl)fe2 OM; A) = D%D%<j>ao(x,t; A) (27) 
and the corresponding moment expansion of a function by 

00 

u(x,t)= M[k 1 ,k 2 ;t](f> klM (x,t;X) (28) 
fei,fc 2 =i 

Note that if the function u(x,t) in (128]) is the vorticity field of some fluid the 
linearity of the Biot-Savart law implies that we can expand the associated velocity 
field as: 

oo 

V(z,t)= ^i^2;t]V fel , fc2 (x,t;A) (29) 

fcl,fc2=l 

where 

V klM (x,t;\) = D k x lD^V m (x,t;X) (30) 
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and V o(x,t; A) is the velocity field associated with the Gaussian vorticity distribu- 
tion 0oo - explicitly we have: 

V o(x,t;A) = ^t^(l-e-^ 2 ), (31) 

We define the Hermite polynomials via their generating function: 

H ni , n2 (z; A) = ((D%D%e(^)) \ t=0 (32) 

Note that the "standard" Hermite polynomials correspond to taking A = 1. 

Then using the standard orthogonality relationship for the Hermite polynomials: 

/ H nun2 (z; A = l)H mi>m2 {z; A = l)e~ z ' 'dz = 7r2 ni+ri2 (n 1 !)(n 2 !)5„ 1 , ml (5 n2 . m2 (33) 
Jm? 

we see that the coefficients in the expansion fl28|) are defined by the projection oper- 
ators: 

M[k u k 2 ;t] = (P klM u)(t) = — ^ j^H kl!k2 (z;X)u;(z,t)dz (34) 
3.1 Convergence of the moment expansion 

In this subsection we derive a criterion for the convergence of the moment expansion 
derived above and we show that if this criterion is satisfied for t = then it is satisfied 
for all subsequent times t > 0. 

Our convergence criterion is based on the observation that the Hermite functions 
0fc lj / C2 (x, t; A) are, for any value of t, the eigenfunctions of the linear operator 

C x ^j = ^A 2 At/; + • (x-0) . (35) 

This fact can be verified by direct computation and is related to the fact that C x 
can be transformed into the Hamiltonian quantum mechanical harmonic oscillator. 

The Gaussian function 0oo plays a crucial role in the convergence proof, and its 
dependence on the parameter A is particularly important in this discussion, so for 
this subsection only we will define 

0ooOM; A) = $\{x,t) 

to empasize this dependence. 

If one now proceeds as in Lemma 4.7 of [5] one can prove 
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Proposition 3.1. The operator £ is self- adjoint in the Hilbert space 

X x = {f e L 2 (M 2 ) | $" 1/2 / G L 2 (M 2 )} 
with innerproduct (f,g)\ = L 2 ^^fgdx. 

An immediate corollary of this proposition and the general theory of self-adjoint 
operators is 

Corollary 3.2. The eigenfunctions of C x form a complete orthogonal set in the 
Hilbert space X x . 

and as a corollary of this result and the observation that the eigenfunctions of C x 
are precisely our Hermite functions {4>ki,k 2 }i we have finally 

Proposition 3.3. Suppose that 

\\f\\l= [ $ A " 1 (x)|/(x)| 2 ^<oo, 

then the expansion 

f(x) = ^M[k u k 2 ]<f> klM {x) 
fci,fc 2 

converges with respect to the norm on the Hilbert space X x . 

Thus, the following criterion guarantees that the expansion (|28|) for the vorticity 
converges: 

/ ^ 1 (x)(u(x,t)) 2 dx < oo . (36) 
Jr 2 

The main result of this subsection is the following theorem which proves that if 
our initial vorticity distribution satisfies fl36l) for some A = Ao, then the solution of 
the vorticity equation with that initial condition will satisfy (136]) for all time t with 
A = a/ 4z4 + Ag and hence as a corollary if the initial vorticity distribution satisfies 
( 1361) then our moment expansion converges for all times t. 

Theorem 3.4. Define the weighted enstrophy function 

£(t) = [ $- x (x)(w(x,t)fdx . 

If the initial vorticity distribution Uq is such that S(0) < oo for some Ao, and uq is 
bounded (in the L°° norm) then S(t) is finite for all times t > 0. 
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Proof: The idea of the proof is to derive a differential inequality for 8 (t) which 
guarantees that if 8(0) is finite then 8(f) will be finite for all t. Differentiating 8(t) 
we obtain 

= %m - Y.jjx^^x^xMdx (37) 
+2/ ^~ 1 (x)u(x,t)d t u(x,t)dx 

JR 2 

Au Au r 

= jpS® - Yi Jjx\^(x)(^t)fdx (38) 

+2 ^ 1 (x)u(x,i) (uAlj - u ■ Vlj) dx 
Jr 2 

We now consider the last term in f)38p . First note that upon integration by parts 
we have: 

2 / ^(x)uj(x,t) (vAou(x,t))dx= -2v [ ^(x) ( \Vu\ 2 + ^-ux ■ Vu) dx . (39) 

JR 2 J R 2 \ A J 

The right hand side of the expression in (1391) can again be broken up into two pieces 
and the second can be bounded by: 

2v I ^(x) f^wi ■ Vw] dx<v\ §^(x)\Vu)fdx + / $^ (xV) dx . (40) 
Jr 2 \a J Jr 2 A J R 2 



Finally, we bound the last term in (1371) . which comes from the nonlinear term in 
the vorticity equation. In this estimate we use the fact (see [2], Lemma 2.1) that 
the L°° norm of the velocity field u can be bounded by a constant times the sum 
of the L l and L°° norms of the vorticity field - i.e by C(||u;||£i(]R2) + HwHioo^)). 
This observation, combined with the fact that £)||lj>(r 2 ) < ||c<jo||lp(ir 2 ), which is 
a consequence of the maximum principle, implies that 

II u ('^)I|l°°(r 2 ) < Cdl^olU 1 ^ 2 ) + ||u;o||l°°(I8 2 )) • 

and hence that 

2/ <5>^(x)uj(x,t) (u- Vu)dx < 2C(u ) / ^ 1 (x)\u(x,t)\\Vuj\dx (41) 
Jr 2 Jr 2 

^4CVq) f $-i( a .)( a; ( a . ji ))2 da . + I/ f §-\ x )\Vuj\ 2 dx 
v Jr 2 Jr 2 
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If we now combine the inequalities in ( J39l) . (1401) and ( 14T1) . with the expression for 
§ in (J3ED we obtain: 

| (i) ,(^ + -) £(t) , (42) 

from which we see immediately that if £{0) is bounded, remains bounded for 
all time. □ 

3.2 Differential equations for the moments 

Assuming that the function u(z,t) is a solution of ([7]), we can derive differential 
equations satisfied by the moments M[ii,k 2 ,t] in ([28]) . Surprisingly the expressions 
for the coefficients in these expansions are quite simple and explicit. 

If we differentiate ( 1281) and assume that u is a solution of the two-dimensional 
vorticity equation we obtain: 

°° cLAdlk k ' t\ °° 
d t uj = ^2 J t 2 ' #fci,fea(M;A)+ ^I^i> h]t]d t (pk u k 2 ( x ^^ ^) 

fel,fe2=l fcl,fc2 = l 

OO 

= ^ M[fci,fc 2 ;t] (z/A0 fel)fe2 (a:,t; A)) (43) 

fci,A;2=l 

(oo \ / oo 

M[e lt £ 2 ;t]Vt l<t3 (x t t;\) I • V I ^ Mf^, A; 2 ; t]0 fel , fc2 (a;, t; A) 
£i,£ 2 =i / \fci,fe 2 =i 

From the first of the "simple facts" we stated about 0oo ; we see that the last term 
on the first line cancels the middle line and hence if we apply the projection operators 
defined in (1341) . we are left with the system of ordinary differential equations for the 
moments 



dM[h,k 2 ;t] 

- — —ii 



dt 



ki,k 2 



J2 4; t\V ilt e 2 (x, t; A) 1 (44) 

(OO 
^ M[m 1 ,m 2 ;t](f) mijm2 (x,t; A) 
mi,m 2 =l 

The somewhat surprising fact which, in our opinion, makes the preceding straight- 
forward calculations interesting is that the projection on the right hand side of (jUJ) 
can be computed explicitly in terms of the derivatives of an relatively simple function. 
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We now explain how this is done. First recall that: 

(f> mi ,m 2 (x, t- A) = D^D™><f> 00 (x, A), V h , h (x, t- A) = D e ^D%V 00 (x, A) (45) 
In order to avoid confusing the two sets of derivatives we will rewrite these formulas 



as 



n>m2 (x,t;X) = (D^DZ 2 cj>o (x + a,\))\ a=0 , (46) 



V* A (:M;A) = (DtlDtlV 00 (x + b,\))\ 



6=0 



Inserting these formulas into the right hand side of (1441) and using the formula for the 
projection operator Pk lt k 2 i n terms of the integration against a Hermite polynomial 
we obtain 

dM ('_]_)(fei+fc2)^2(fci+fc 2 ) 

— [A^WH- ^ ; J2 E M[£ u £ 2 ,t]M[ mil m 2l t] 

x f H klM (x)(D^D^V 00 (x; A)) ■ V x (D^> (x; A))dx 

(_l){kl+k 2 )\2(ki+k 2 ) 

= ~ ' 2 L 2{kll)m E E M[^ 2 ,t]MK,m 2 ,t] 



f-l,t2 mi ,"12 

* D ti D % D Ti D T 2 D ai D a 2 ^a- [ I e ' ' ^' V 00 (x + 6; A)0 O o (* + a; A)<fa 

The last equality in this expression results from rewriting the Hermite polynomial 
Hi lt k 2 in terms of its generating function. 

The last step in deriving the equations for the moments is to evaluate the integral 
in the last line of (1471) . The key step in this evaluation is to recall that for these 
incompressible flows the velocity field can be written in terms of the derivatives of the 
stream function and that the Laplacian of the stream function is minus the vorticity. 
Thus, we can write: 

V 00 (x + b;X) = -V£(A 6 )-Voo(z + 6) , (48) 

where V*J = (d x J,-d x J). 

Inserting this into the integral in (147|) we find 

(-t 2 ^-t 2 .+ 2 t 1 x 1 + 2 t 2 x 2 ) 

e & V 00 (x + 6; \)(p 00 (x + a; \)dx (49) 

r (— ( i— t^+ 2 t 1 x 1 + 2 t 2 x 2 ) 
= -V* b (A b )- 1 / e 0oo(x + M)0 oo (x + a;A)cix 

Jr 2 
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(-t{-tl+ 2 t 1 x 1 +2t 2 x 2 ) 



Now note that all three factors in the integrand are Gaussians and thus the integral 
can be evaluated explicitly, and we find 



e a 2 </>oo(^ + b; X)4>oo(x + a; X)dx. (50) 

_ __}_ e ~^( K al+a^~2a 1 b 1 +b 2 L ~2a 2 b2+b^+2ait 1 +2b 1 t 1 +tj+2a2t2+2b2t2+t^) 



2ix\ 2 

We next compute the expression 

V*(Ab)~ 1 - c -^( a ?+ a ^- 2a i fe l+ b 1 ! - 2a 2 fe 2+^+2aiti+2&ifi+^+2a2t2+2b2t2+^) 
27rA 2 

Recall that given a vorticity field to, — (A)~ l io is the associated stream function and 
— V(A) _1 u; the velocity field associated with to. Since the inverse Laplacian and 
derivatives in floTT) act only on the 6-dependent parts of the expression we need to 
evaluate 

V*(Afe)~ 1 c -^r(-2ai6i+6?~2a 2 62+6j*+2Ml+2fe2*2) 

b 2n\ 2 

= _ e5 ^((il-«i) 2 +(t2-a 2 ) 2 )y*^^-l^_ e - 3 ^((fe 1 +(ti-ai)) 2 + (fe 2 +(t 2 -a 2 )) 2 ) 

But 

_y*^ b )-l J_ e -^2-((6l + (*l-«l)) 2 + (fe 2 + (* 2 -a 2 )) 2 ) 

is just the velocity field associated with a Gaussian vorticity distribution (i.e. an 
Oseen vortex) centered at the point — ((t\ — ai), (t 2 — 02)) which we know explicitly. 
Hence, the expression on the right hand side of (1491) has the explicit representation: 

J_ e -5^(«?+»i+*?+*i) e ^((*i-»i) 2 +(fe-a2) 2 ) ( 53 ) 

(-(62 + (t2-a 2 ),(6l + (tl-ai)) / , ^(( 6l +(t 1 -a 1 ))a + (6 2 +(t 2 -a 2 ))a)\ 

((b 1 + (t 1 -a 1 )) 2 + (b 2 + (t 2 -a 2 )) 2 ) I 7 

If we now return to (|47|) we see that in order to compute the coefficients in the 
moment equations we need to evaluate the divergence of this last expression with 
respect to a which gives 

K(a 1 ,a 2 ,b 1 ,b 2 ,t 1 ,t 2 ;X) = 

(a 2 h -b 2 ti + (-ai + fei) t 2 ) 
n (a 1 2 + a 2 2 + 6 1 2 + 6 2 2 + 2 6 1 ti+t 1 2 -2ai (61 + h) + 2 b 2 1 2 + t 2 2 - 2 a 2 (b 2 + t 2 )) X 2 

(-a 1 +bl+t 1 ) 2 + (-a 2 +i) 2 + t 2 ) 2 \ ai 2 + a 2 2 -2ai 6 1 +6 1 2 - 2 a 2 i 2 + f, 2 2 +2 ai t 1 +2b 1 t 1 +t 1 2 +2a 2 t 2 + 2 b 2 t 2 +t 2 2 

-1 + e e ^ (54) 
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Returning to equation (l47j) we finally conclude that 



dM 




(55) 



dt 



2 k ^(k 1 \)(k 2 \) 



^[h,k2]£xj2,m 1 ,m 2 -,X}M[£ 1 ,£ 2 ,t}M[m 1 ,m 2 ,t] 



£1/2 m\,m2 



where 



T[k 1 ,k 2 ;£ 1 ,£ 2 ,m 1 ,m 2 ] X] = 

= D^D^D^D^Di\D^K( ai , a 2 , h, b 2 , t u t 2 ; A)| 



t=0,a=0,b=0 



(56) 



Thus, we have succeeded in rewriting the two-dimensional vorticity equation as a sys- 
tem of ordinary differential equations with simple, quadratic nonlinear terms whose 
coefficients can be evaluated in terms of derivatives of a single explicit function. Fur- 
thermore, we have given a sufficient condition on the initial vorticity distribution 
to guarantee that the expansion of the vorticity generated by the solution of these 
ordinary differential equations converges for all time. 

4 The moment expansion for several vortex cen- 
ters 

In this section we extend the Hermite moment expansion of the previous section to the 
case in which there are two or more centers of vorticity by combining this expansion 
with the multi-vortex representation of Section [2J For simplicity of exposition we 
limit the discussion here to the case of two vortices but the expansion can be extended 
to any finite number of vortices. 

The basic idea is just to consider the equations f|T5|) for the evolution of each 
vortex and then expand each of the functions uji in Hermite moments as in the 
previous section. Thus, we define 



for j = 1,2. We make a similar expansion for the velocity field in terms of the 
functions V^^ 2 , and insert the expansions into f[T5l) . Letting z = x — x^(t), and 



00 




(57) 



fcl,/C2=l 
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recalling that dt^M = ^^■4>k 1 ,k 2 we obtain: 

dM j [k 1 ,k 2 ;t] 
dt ~~ 



(58) 



-Pk lt k 2 



2 oo 

(oo 
M j [m 1 ,m 2 ;t}<j) (z,t;X) 
k lt k 2 =l 

where Sjj/ = x J '(t) — x^(t). Proceeding as in the previous section and using equation 
(1471) we finally conclude that 



dM j 
dt 



[fa, k 2 ,t] 



(59) 



.]_^(fci+fc 2 )^2(fc 1 +fc 2 ) ^ 

Z2z2 a2 rJ ' 3 '[fa^ 2 ]hJ2,m 1 ,m 2 ; s jtf , X] 



2 k i+ k *(k 1 \)(k 2 \] 



j'=l m\,m,2 



where 



r J ' 3 '[fa, k 2 ;l 1 ,£ 2 ,m 1 ,m 2 ; s jd >, A] = 

= D k t l D% D™ 1 D™ 2 D £ a \ D% K multi (a u a 2 , b u b 2 , t u t 2 ; s hf , A)| t=0 , a =o,b=o 



(60) 



and 



K multi (a 1 ,a 2 ,b 1 ,b 2 ,t 1 ,t 2 ;s 1 ,s 2 ,X) = (61) 
V a • ( / e a 2 V 00 (x - s + b; X)(poo(x + a; A)dx 



Remark 4.1. Note that comparing ( 160]) with the calculation leading ( ]54|) we see that 
K multl can be written in terms of the expression for K via the simple formula: 



K mu l (a x , a 2 , b x , b 2 , h, t 2 ; s u s 2 , A) = K(a u a 2 , h - si, b 2 - s 2 , t u t 2 ; A) 

We now derive a similar expansion for the evolution of the centers of each vortex. 
We begin with 



dxi 
~dt 



N 

(t) = — V / (u e (z + x J (t)-x j '(t),t)uj j (z,t))dz 
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(Recall that this is really a pair of equations, one for each component of xK) Now 
insert the moment expansion of u J and uji into this expression and we obtain: 

dx j N f 

"dTW = Z)Zl Yl M'[£ 1 ,e 2 ,t]M 3 '[m 1 ,m 2 ,t] / V il;h (z - s j;j ,, t)<f> kl;k2 (z, t)dz , 

(62) 

where as before sjji = x J (t) — x^(t). The integral of the velocity and vorticity can 
be evaluated just as in the preceding section and we find 

d j N 

-^-{t) = J2J2J2 E^\£ 1 ,£ 2 ,m l ,m 2 -,s jJf ,X}M j [£ 1 ,£ 2 ,t}M''[m 1 ,m 2 ,t} , (63) 

j'=l l\,l-x mx,m% 

where in this case the coefficients B 7 '- 3 ' [k±, k 2 ; £x, £ 2 , mi, m 2 ; Sjj>, A] is given by the 
expression 

S^'[4,4,m 1 ,m 2 ; S ,A] (64) 

- J_ r>™i n™2 n^i rfr ( ( ~ s 2 ~ a 2 ), (6i - Si - a - 1)) x 
" 27r ^ ^ ^ (fea _ ^ _ fl2)2 + (6i _ si _ ffl _ 1)2 J 



1 — e 2A 



1 ((b 2 _ S2 _ a2 )2 + (bl _ Sl „ a _ 1 )2 ) 



a=0,b=0 



4.1 Convergence of Multi- Vortex Expansion 

We note that it is easy to extend our previous result on the convergence of the 
moment expansion to the multi-vortex expansion. To do so, we change variables to 
re-center each vortex at the origin. Thus if we let 

J (x - x j (t) , t) = w j (x,t), u j (x - x j (t) , t) = v j (x, t) 

then equation (fT51) becomes 



dw j 



dt 



,r. 



t) = vAw j (x, t) - ^ ( x ~ + xj (t)^) ■ Vw j (x, t), 



=i 



forj = l,...,iV. 

We are now ready to state our result. 
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Theorem 4.2. Define 



£3[t) = / G x (x)w 3 (x,t) dx. 
Jv 2 

If the initial vorticity distribution uo = Y^e=i ^(x — x e (0), 0) is such that £3(0) < oo 
for some Xq, and for all j = 1,2, ... ,N, and if uj$ is bounded (in the L°° norm) then 
each £ J (t) is finite for all times t > 0. 

Proof. We use the same idea as the single vortex case and differentiate £ J (t): 

d -^(t) = p*(t) - ^ \x\*G- x \x) W \xYdx 

+2 / Gl 1 (x)w i (x,t)(d t w i (x,t)dx 
Jr 2 

Au Au r 

= T2 £< (*)-TI / \A 2 Gl\x)(w\xfdx 

A * JR 2 

+2 f G^w^uAw'- ( Vv £ (x-x £ (t) +x j (t),t) ) -Vw^dx. 
Jr 2 \<=i / 

From here the proof proceeds identically to the proof of theorem 13.31 except we 
must examine the nonlinear term that comes from the vorticity equation a bit 
closer. In general, it is unknown whether w J (-,i) satisfies a maximum principle, 

but since f^)^ v e (x — x e (t) + x ] (t),t)j = u, the solution to (TjQ), it does satisfy a 
maximum principle. Hence, as in Subsection 13.11 we can bound the L°° norm of 
Ylf=i v^(x — x e (t) + x 3 (t),t) by a constant depending only on the initial vorticity 
distribution. Hence we proceed to bound the integral 

/ GjV(|Vv f (i-x'(i)+i J (i),i) J • Vw l )dx <2C(lu ) [ G-VilWldx 
Jr 2 ) J R2 

and thus the rest of the bounds are the same. Again, putting everything together 
we arrive at: 



ir w n* + * r (<) (65) 

□ 
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4.2 Interaction of Gaussian Vortices 

The experimental and numerical work of [9] has shown that widely separated regions 
of vorticity can be well approximated by Gaussians for long periods of time. This 
corresponds to truncating our expansions for the vorticity so that they contain only 
a single term. In this subsection we analyze the equations that result from this 
truncation. This approximation can be viewed as a generalization of the Helmholtz- 
Kirchhoff approximation in which we include the effects of vorticity and finite core 
size to lowest order. We show that the total vorticity of each of the vortices is 
constant while the centers of vorticity evolve either along straight lines or circles. 

Remark 4.3. It should be noted here that the effect of truncating our expansion after 
one term allows for only viscosity as the driving force in vortex merger. Allowing 
more terms in the expansion introduces convective forces. 

If we start with two Gaussian distributions for our initial vorticity with the same 
value of Ao , and truncate the equations of motion for the moments so that all terms 
containing higher order moments are omitted then we can conclude that the two 
vortices travel along circular or straight line orbits around the "center of vorticity." 
Moreover, the leading coefficients M J [0,0 : t] of the expansions are constant. To be 
precise, we let 

u\x,t) = M 1 [O,O;t]0 o ,o(a?-a? 1 (*),*;A) (66) 
u 2 {x,t) = M 2 [0,0;t}<f)o fi (x - x 2 {t),t;\). (67) 

Let us also write s i)3 - = x^(t) — x l (t), then using (1551) - (l6Tj) we first calculate the 
evolution of 1^(0,0 : t) 

—M 1 [0, 0; t} 2 K(0, 0, 0, 0, 0, 0, A) 

-M^O, 0; t]M 2 [0, 0; t]K multi (0, 0, 0, 0, 0, 0, s 1>3 , A) 
0, 

since K(0, 0, 0, 0, 0, 0, A) = and K multi (0, 0, 0, 0, 0, 0, A) = 0, respectively. The 
calculation for M 2 [0,0;t] is the same. Thus the evolutions of the coefficients of 
leading order are constant. 

Remark 4.4. The fact that M l [0, 0; t] is constant in time is not a consequence of the 
truncation of the moment equations to first order. One can show that the equations 
(Tl4"l) conserve the zeroth moment of uji , independent of any truncation. 



dM 1 
dt 



[0,0; t] 



19 



More interestingly, though, is the calculation for the evolution of x J '(t), the centers 
of these vortices. Again if we denote x^ = (x?i,2? 2 ) then the evolution for each 
component as defined by equation (1231) can be written as: 

dxi i 



dt M* [0,0, t] 



J Vf(y-s 1>2 ,t)u j (y,t)dy. (68) 



Now using equation (1481) to evaluate 

V k (j — si t 2,t), equation (1681) yields the following equations for each x\, j — 1,2 and 
z = 1,2: 

(x^-^f) 2 + (4-3:|) 2 

•1 M 2 (e 2A (^ 

X l 2tt (xj-x 2 ) 2 + (xi-x 2 ) 2 

(a;-^ — x ) -|- (a^ — ^ 2 ^ 

•1 M 2 (e~ 2A <') 2 -l)(-xi+x 2 ) 

X 2 27T ( X i_^)2 + (x i_^ ) 2 

(,j-.?) a +( ! 4— g) a ( 69 ) 

2 A/ 1 (g 2A ^ -l)(^-s§) 

X l 2tt (x}-x 2 ) 2 + (i;i-a; 2 ) 2 

/ 1 2^2 I / 1 2 \ 2 

•2 A/ 1 (e~ 2A M" -l)(-xl+x 2 ) 

X 2 2tt (x}-a; 2 ) 2 + (a:i-x 2 ) 2 

where M J = M J [0, 0, i] is constant and represents the total vorticity of the j th vortex. 

Remark 4.5. If the vortices have different X(t) values say X\(t) and \2{t) then one 
just needs to replace the 2X(t) 2 with Xi(t) 2 + X2(t) 2 in the exponential to arrive at the 
correct system. 

We now state the main result of this section. 

Theorem 4.6. System ( fgPj) admits only circular or straight line trajectories. 

Proof: Away from rest points our system f)69p can be transformed to: 



dxi 
dx 

dx\ _ m 



2 1 2 2 

2 —^i+^i 

„1 „,2 



I 3 (70) 

dxj M 1 

dx\ _m 2 

dx\ M 1 ■ 
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Integrating out the bottom two equations we get 

x i = — a7t( x i + ^1) 



(71) 



First let us assume that M 2 7^ — M . Then plugging back into the first two 
equations of (1701) we arrive at: 



34 _ (-i-frH-fc 



9x\ (i+M^^i.^ x\-k 2 



(72) 



for appropriate constants k 1 , k 2 , k 1 , k 2 - Thus we integrate again and get 

{x l^_~ k2) 2 +{x l {t) _~ kl) 2 = Ci 

(xl(t) - k) 2 + (xKt) - h) 2 = c 2 . 1 J 

If M 2 = —M 1 then we have equal but opposite size vortices and equations (169]) 
become: 

9x\ _k2 

aS hi / 74 n 

4 — k2 ^ ) 



which gives us straight line solutions with slope — 1|, as desired. 



dx\ _k2 
dxf kl 



□ 

2 



If we now consider the case of equal total vorticity for the two vortices (M 
M 1 = M) , the classical point vortex result is that the vortices will rotate around the 
center of vorticity at a constant frequency, 2 J£ 2 , where D is the distance between 
the vortex centers. We will now compute the viscous and finite core size effects on 
the frequency of rotation, predicted by our model. For simplicity we will center the 
vortices at the origin and place them on the circle of radius, r. We apply the polar 
change of variables, 

x\ = rcos(6>j) x 2 = rsin(#j) i = 1,2 (75) 

then using the fact that our vortices are out of phase by it we can compute that 

(x\ — x\) = r cos(#i) — r cos(6* 2 ) = 2r cos(#i) (76) 
(x\ — x\) = r sin(6 l i) — r sm(9 2 ) = 2r sin(#i), (77) 
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and thus we arrive at the expression for the frequency of rotation 

M ~ 2r2 
fi=_(l- e W). (78) 

Since \(t) 2 = X^ + Aut we notice two things: the first is that viscosity, u, slows the 
frequency of rotation down and, second, in the formal limit as v — > and Aq — > we 
recover the constant frequency, 2 ^ 2 , which is what the Helmholtz-Kirchoff model 
for the rotation of two point vortices predicts. 

With this calculation we may now compare the results of expanding a two gaus- 
sian initial distribution as a single vortex or as two independent vortices. In [10] 
a two Gaussian initial distribution a distance 2r apart with core size Ao, each with 
mass M, is approximated by a single vortex expansion using (|2~8"1) . The authors trun- 
cate the expansion to quadrapole moment (n = 2) and calculated the frequency of 
rotation to be: 



8tt 



2ut V r<2 J r^l + Avt/r 2 



(79) 



This equation is directly comparable to (l7"8l) . Notice that both equations for fre- 
quency of rotation indicate slowing of rotation over time, albeit at different rates. In 
addition, both equations recover the Helmholtz-Kirchhoff approximation of 2 ^ 2 in 
the limit as v — > and Ao — > 0. In fact, in the sufficiently localized regime, r << Ao 
both equations have similar initial frequencies and remain close asymptotically. A 
typical example is shown in figure 14.21 



5 Conclusions: 

In this paper we have derived a system of ordinary differential equations whose solu- 
tions give a representation of solutions of the two dimensional vorticity equation in 
terms of a system of interacting vortices. We have also derived a sufficient condition 
on the initial vorticity distribution which guarantees that this representation in terms 
of interacting vortices is equivalent to the original solution of the two-dimensional 
vorticity equation. This model generalizes the classical Helholtz-Kirchhoff model of 
interacting, inviscid, point vortices to include the effects of both finite core size and 
viscosity. We have also looked at the analytical predictions of our model for the 
interaction of two vortices which the expansion is truncated at leading order. We 
plan in future work to further explore the analytical and numerical predictions of 
this model. 
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"Typical" frequency of rotation for two localized vortices 

0.07- 
0.06- 
0.05- 

c 

J 0.04- 
0.03- 
0.02- 
0.01- 

100 200 300 400 500 

time 

Figure 1: Here we plot both the frequency of rotation predicted by single vortex 
expansion up to quadrapole order (dashed) and the frequency predicted by two vortex 
expansions truncated to leading order (solid). The parameter values used are v = 
.01, M = l,r = 1, and A = .01. 
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